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"*^ ; ABSTRACT 

f^ ■ The presence of nearby discrete cosmic-ray (CR) sources can lead to many interesting effects 

CO \ on the observed properties of CRs. In this paper, we study about the possible effects on the 

CR primary and secondary spectra and also the subsequent effects on the CR secondary-to- 
(— I \ primary ratios. For the study, we assume that CRs undergo diffusive propagation in the Galaxy 

pi !■ and we neglect the effect of convection, energy losses and reacceleration. In our model, we 

assume that there exists a uniform and continuous distribution of CR sources in the Galaxy 
\^ • generating a stationary CR background at the Earth. In addition, we also consider the exis- 

jj^ \ tence of some nearby sources which inject CRs in a discrete space-time model. Assuming 

C^ ■ a constant CR source power throughout the Galaxy, our study has found that the presence 

of nearby supernova remnants (SNRs) produces noticeable variations in the primary fluxes 

mainly above ~ 100 GeV/n, if CRs are assumed to be released instantaneously after the su- 
^ I pernova explosion. The variation reaches a value of ~ 45% at around 10^ GeV/n. Respect to 

earlier studies, the variation in the case of the secondaries is found to be almost negligible. We 

also discuss about the possible effects of the different particle release times from the SNRs. 
^-1- i For the particle release time of ^ 10^ yr, predicted by the diffusive shock acceleration theories 

in SNRs, we have found that the presence of the nearby SNRs hardly produces any significant 

effects on the CRs at the Earth. 
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^ ' 1 INTRODUCTION 

H , 

^ Cosmic-rays (CRs) below the knee region (~ 3 x 10^^ eV) are generally considered to be of Galactic origin. Though the nature of their 
sources are not exactly known, it is widely beleived that majority of them are accelerated in supernova remnant (SNR) shock waves. Such 
a hypothesis is mainly based on the similarity of the power supplied by a supernova explosion to the power required to maintain the CR 
energy density in the Galaxy (~ 10''^ erg s^^). In addition, studies using nonlinear effects produced by accelerated CRs have shown that the 
maximum energy of CRs that can be accelerated in SNRs is ^ Z x 10^^ eV (where Z is the total charge of the particle) which is close to 
the knee (Berezhko 1996). 

It is quite expected that the presence of one or more nearby discrete sources can result in various significant effects on the observed 
properties of CRs. Studies concerning their contributions to the total CR electron flux can be found in earlier studies like Mao & Shen (1972) 
and Cowsik & Lee (1979) under the framework of an energy independent CR diffusive propagation model and later, in Atoyan et al. 1995 
and Kobayashi et al. 2004 assuming a more realistic energy dependent diffusion cofficients in the Galaxy. For CR nuclei, Lingenfelter (1969) 
studied the density variations at the Earth due to the local sources using a simple particle escape model and later, Erlykin & Wolfendale (2001) 
carried out Monte Carlo calculations assuming SNRs to be distributed stochastically throughout the nearby Galaxy. Strong & Moskalenko 
(2001) studied the influence of the discrete nature of SNRs on the CR proton densities in the Galaxy using their GALPROP CR propagation 
code. In addition, some recent studies tried to stretch the importance of local sources to explain some strongly observed features of CRs like 
the knee and the anisotropy. For instance, Erlykin & Wolfendale (2000) claimed that the knee in the CR spectrum can be attributed to the 
presence of a single recent supernova (as yet unidentified) in the local region. Erlykin & Wolfendale (2006) further tried to explain the rise 
in the anisotropy amplitude as well as the change in its phase near the knee using a single source exploded in the direction from the Sun 
downward of the main CR flux, which comes predominantly from the inner Galaxy. Moreover, using a convection-diffusion CR propagation 

* E-mail: satyend@barc.gov.in 



2 S. Thoudam 

model, Thoudam (2006) tried to emphasise the necessity for the presence of at least one old nearby source in order to explain the observed 
proton flux below ~ 100 GeV/n. In another recent study, Thoudam (2007a) studied the effect of the nearby known SNRs on the observed 
CR anisotropy below the knee and later on, suggested the possibility of explaining the data above 100 GeV/n by a single dominant source 
with properly chosen source parameters and claimed that the source may be an undetected old SNR with a characteristic age of --^ 1.5 x 10^ 
yr located at a distance of ~ 0.57 kpc from the Earth (Thoudam 2007b). 

Taillet & Maurin (2003) pointed out that the majority of the CRs reaching the Earth are possibly emitted by sources located within a 
few kpc and hence, the propagation parameters that we usually derived from the observed secondary/primary (s/p) may only give a local 
information which may be different from other parts of the Galaxy. In another work, Biisching et al. (2005) showed that even if the CR 
propagation parameters are assumed to be constant throughout the Galaxy, the discrete nature of the CR sources can produce a much larger 
fluctuations of the CR primary densities than that of the secondaries, implying a significant fluctuations in the s/p ratio. In this paper, we shall 
also study about the expected fluctuations on the CR primary and secondary spectra at the Earth, but by incorporating the nearby known CR 
sources (considered here as SNRs). Such an approach is expected to give a more detailed understanding about the CR density fluctuations 
and its implications as far as the position of the Earth is concerned. 

We plan the paper as follows. In section 2, we give a brief description about our model. In section 3, we calculate the CR primary 
and secondary spectra from a discrete point-like source and in section 4, we calculate for the CRs from a continuous and stationary source 
distribution. In section 5, we study the CR density variations at the Earth due to a single nearby source and in section 6, we give an 
application to the nearby SNRs and present a comparison of the calculated CR spectra with the observed data. Finally in section 7, we give a 
brief discussion about our results and its implications. 



2 MODEL DESCRIPTION 

Our model assumes that a major fraction of the CRs detected at the Earth are liberated from sources which are distributed uniformly and 
continuously (both in space and time) in the Galaxy. We refer to these sources as the background sources and the CRs they emit as the 
background CRs. In addition, we also assume that there exists some nearby discrete sources which inject CRs in a discrete space-time model 
and whose contributions to the total CRs are yet to be investigated. The contributions of the discrete and the background sources will be 
treated seperately. CRs from the discrete sources will be discussed under the framework of a time dependent diffusive propagation model 
while those from the background sources in a steady state model. 

For the background CRs, the flattened shape of our Galaxy allows us to assume their diffusion region as a cylindrical disk of infinite 
radius with finite half-thickness H. Such an assumption of infinite radius is valid at least for CR studies at the position of the Earth, which is 
at a distance of ~ 8.5 kpc from the Galactic Centre, where the effect of the Galactic radial boundary (assumed to be > 20 kpc) on the CR 
flux is expected to be negligible. This is because a substantial fraction of the CRs reaching the Earth are liberated from sources located within 
an approximate distance which is of the order of H (see e.g. Taillet & Maurin 2003). The actual value of H is not exactly known. The values 
estimated using different CR propagation models fall in the wide range of (2 — 12) kpc (see e.g., Lukasiak et al. 1994, Webber & Soutoul 
1998, Strong & Moskalenko 1998). For our calculations, we choose a value of H — ±5 kpc. We further assume that the background sources 
as well as the interstellar matter are distributed in a thin disk of radius R and half-thickness h. Here again, the effect of R on the CR flux is 
negligible as long as i? > H. Observations have found that both the distributions of the SNRs as well as those of the atomic and molecular 
hydrogen extend approximately upto a radial distance of R ^ 16 kpc (Case & Bhattcharya 1996, Gordon & Burton 1976). Regarding the 
vertical distributions, detailed studies have found that most of the SNRs and the molecular hydrogen are confined within the region ~ ±200 
pc from the Galactic plane (see e.g Stupar et al. 2007, Bronfman et al. 1988). The distribution of atomic hydrogen also follow a similar 
structure but with a thin long tail extending as far as ^ 700 pc from the plane (Dickey & Lockman 1990). Therefore, in our calculations we 
consider an infinitely thin disk approximation for both the background sources and the matter distributions. Such an approximation can also 
be found in some earlier works like Webber et al. 1992, Seo & Ptuskin 1994 etc. 

For the nearby discrete sources, we consider only known SNRs located within a distance of 1.5 kpc from the Earth. They are listed in 
Table 1 along with their estimated distances and ages. They are the sources which are expected to produce significant temporal fluctuations 
in the CR densities at the Earth (Thoudam 2006). The diffusion region for CRs from these sources is assumed to be of infinite dimensions. 
This assumption is based on the fact that CRs from nearby sources do not effectively see the presence of any Galactic boundaries (both in the 
radial as well as vertical directions) due to their much smaller propagation times to the Earth compared to their escape timescales from the 
Galaxy boundaries (Thoudam 2007b). 

One more assumption that we make in our model is that since the Earth is reported to be only ~ 15 pc away from the Galactic median 
plane (Cohen 1995), we will simply consider that the Earth is located on the median plane itself in all our calculations. 
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Table 1. Parameters of known SNRs located within a distance of 1.5 kpc from tlie Eartli (See tfie references given in Thoudam 20()7a). 



SNR 


Distance(A;pc) 


Age(yr) 


G65.3+5.7 


1.0 


14000 


G73.9+().9 


1.3 


10000 


Cygnus Loop 


0.4 


14000 


HB21 


0.8 


19000 


Gl 14.3+0.3 


0.7 


41000 


CTAl 


1.4 


24500 


HB9 


1.0 


7700 


S147 


0.8 


4600 


Vela 


0.3 


11000 


G299.2-2.9 


0.5 


5000 


SN185 


0.95 


1800 


Monogem 


0.3 


86000 


Geminga 


0.15 


340000 



3 CRS FROM A DISCRETE POINT SOURCE 
3.1 CR primaries 

In the diffusion model, neglecting convection, energy losses and particle reacceleration processes, the propagation of CR primaries in the 
Galaxy can be represented by the equation 

8N 
V • (DpVNp) - 2hnVpapS{z)Np + Q^ = __p (1) 

where the subscript p denotes the primary nuclei, Np{r, E, t) is the differential number density at a distance r at time t, E is the kinetic 
energy per nucleon of the nuclei, Dp{E) is the diffusion coefficient and n is the target density both of which are assumed to be constant 
in the Galaxy, ap is the primary spallation crossection assumed to be independent of energy, Vp is the primary velocity and Qp(r, E, t) — 
Qp{E)S{r)S{t — to) is the particle production rate from the source. In Eq. (1), we neglect the yield of the primaries from the fragmentation 
of heavier nuclei. 

We are also interested in the study of the secondary nuclei which are produced by the nuclear interactions of the primaries with the 
interstellar medium (ISM). Since the energy per nucleon is an almost conserved quantity in the fragmentation process, we will be dealing 
with the energy per nucleon in all our relations rather than the total kinetic energy of the nuclei. Using Green's function technique and 
performing proper fourier and laplace transforms, the exact unbounded solution of Eq. (1) for a source having particle release time to is 
obtained as (see Appendix A) 

r ^2 1 

87rZ3^/'(t-to) 1 VAi - *o V 2^/r^J < 

where q — \z\/y^ Dp , b = 2hncap/2y^Dp and we have taken Vp ~ c, the velocity of light since we are dealing with high energy particles. 
The particle flux can then be calculated using /p(r, E, t) ~ (c/47r)A'p(r, E, t). It should be noted that we assume the source spectrum to be 

Qp{E) = Apqp{T) with qp{T) given by 

qp{T) = k{T^ + 2Tmp)-'^''+^^''^{T + nip) (3) 

where T — ApE represents the total kinetic energy of the nuclei, Ap is the mass number, mp is the mass energy, F is the spectral index and 
k is the normalization constant. 



3.2 CR secondaries 

The transport of CR secondaries in the Galaxy also follow an equation similar to that of the primaries given above as 

V • [DsVNs) - 2hnvsOs5{z)Ns + Q- = ^ (4) 

where the subscript s represents the secondary nuclei and all the quantities have the similar definitions as in Eq. (1). In our model, we assume 
that the secondaries are the results of fragmentation of one or more heavier primaries and we neglect the production of secondaries at the 
source. For secondaries originating from a single type of primary, we can write the source term in Eq. (4) as 

roc J 

qs{r,E,t) = 2hnc5{z) ■j^aps{E,E')Np{r,E' ,t)dE' (5) 
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We can approximate the differential production crossection daps {E, E')/dE' of an s-type nuclei of energy per nucleon E by the fragmen- 
tation of a p-type nuclei of energy per nucleon E' by a delta function as 

-^'^PsiE, E') = apsS{E' - E) (6) 

where Ops denotes the total fragmentation crossection of p to s. This simplifies Eq. (4) as 

Qs(r, E, t) = 2hnc(jps5{z)Np(r, E, t) (7) 

Note that here Np{r, E, i) is given by Eq. (2). For a secondary source term given by Eq. (7) we can easily obtain the solution of Eq. (4) at 
the spatial location r = as 



Ns (E, t) — 2hncaps I dto j dl^cxp 



(r-r')^ 



ADsit-to) 



2hncas 



87r3/2[75,(t-to)]'/' 16Di(i-io) 



exp 



4D, 



xerfc 



2hncas\/t — to 



2Vd: 



Np{r',E,t)S{z') (8) 



where r' denotes the position of the primaries with respect to the point source which in cylindrical coordinates yield an integral over the 
volume element as 

/■oo p27T roci 

dV ^ r'dr' / d<^' / dz' (9) 

Jo Jo J-oo 

Eq. (8) gives the CR secondary density at r = 0, which is considered here as the position of the observer, due to a single type of primary 
emitted from a point source located at a distance r from the observer. If a particular secondary species is produced as the result of fragmen- 
tations of more than one type of heavier primaries, the total number of secondaries produced is obtained by simply adding the contributions 
from the different primary species as 

nI'"^J2^s (10) 

where j denotes the primary species. 

4 CRS GENERATED BY THE BACKGROUND SOURCES 

4.1 CR primaries 

CRs from the background sources are assumed to follow a steady state propagation equation in the Galaxy as given below 

V ■{DpVNp)-2hnVpap5{z)Np^ -Sp (11) 

where Sp{r, E) represents the source term. For the presence of a Galactic vertical boundary at z — ±H, the solution of Eq. (11) for a 
uniform source distribution Sp{r, E) = '^Qp{E)5{z) extended upto a radial distance R in the Galactic plane is given by (see Appendix B) 

,, , ^, R'StQp(E) f^ sinh[K(H-z)] ■r,r^r.^,r. 

Np{z,E)= ^^^ ' / p^^ L^ jXj,{KR)dK (12) 

2-Dp Jo smh{KH) \Kcoth{KH) + ^^^ 

where Ji is the Bessel function of order 1 and 5R = 25 Myr^^ kpc^'^ denotes the supernova explosion rate in the Galaxy (Grenier 2000). 
As already mentioned in section 2, a major fraction of the CRs reaching the Earth are liberated by sources located within a distance which 
is of the order of the vertical height H — ±5 kpc (see Taillet & Maurin 2003). Therefore, considering the fact that the Earth is located at a 
distance of ^ 8.5 kpc away from the center of the Galaxy and that the sources are distributed upto a radial distance of i? --^ 16 kpc, Eq. (12) 
can be used to obtain the CR density at the Earth by setting z = 0. 

4.2 CR secondaries 

The density of CR secondaries which are produced by the fragmentation of the heavier primaries originated from a stationary source distri- 
bution can be obtained similar to that of the primaries as 

Ns{z,E) = 2hncapsNpiz,E)^ H smh[K{H - z)] x Jj{KR)dK (13) 

^^' Jo sinh{KH) Kcoth{KH) + ^^ 

where we have taken the secondary source term as Ss{r, E) — 2hncaps5{z)Np{r, E). Taking z = Q, Eq. (13) gives the secondary densities 
at the Earth. Eq. (13) shows that for very large Ds, the secondary to primary ratio for a stationary source distribution follows 

The above relation shows that if CRs are liberated by stationary sources which are distributed uniformly, the s/p can be used to estimate the 
CR diffusion cofficient in the Galaxy. However, this approach may fail in cases where the influence of nearby discrete sources is significant. 
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Figure 1. B/C ratio at the Earth due to the background CRs (sohd Hne). Model parameters: Oxygen/carbon source abundance ratio O/C =1.4 and <I> = 500 
MV. Experimental data are taken from Swordy et al. 1990, Panov et al. 2007 and the compilations of different experiments given in Stephens & Streitmatter 
1998. 

It is because the presence of strong nearby sources can significantly affect the CR fluxes (mainly the primaries) which can subsequently affect 
the s/p ratio (see Biisching et al. 2005). How far the measured ratio will deviate from Eq. (14) actually depends on the ages and distances 
of the nearby sources. In the next section, we will investigate the variations that one can expect in the CR densities due to the presence of a 
single nearby source. 



5 CR DENSITY VARIATIONS DUE TO A SINGLE NEARBY SOURCE 

One very important parameter in the study of CR propagation in the Galaxy is the CR diffusion coefficient D. The diffusion of CRs is 
generally considered to be due to scattering either by magnetic field iregularities or by self excited Alfven and hydromagnetic waves. 
Because of the possible spatial variations in the scattering processes, the value of D may be different at different locations in the Galaxy. 
However, in the present work we make a simple approximation that D remains constant throughout the Galaxy. As already discussed in the 
last section, the diffusion coefficient can be determined using the s/p ratio sunder the steady state model. Though the value thus obtained 
may not represent the true value because of the presence of nearby discrete sources, to begin with, we assume that 



A(S») = A) 



Do[^ 



{E, < Eo) 



{E, > Eo) 



(15) 



where Ei denotes the kinetic energy per nucleon of a particular nuclear species denoted by the subscript i and the values of Do and Eo are 
chosen so that Eq. (14) fits the observed boron/carbon (B/C) ratio at 1 GeV/n (see Figure 1). We obtain Do = 2.9 x 10^** cm'^s^^ and 
the particle rigidity (corresponding to energy Eo) with charge Zi and mass number Ai as po — AiE^/Zi — 3 GV. We assume that the 
boron secondaries (^^B) are produced by the nuclear interactions of the ^^C and ^^O progenitors with the ISM. The values for the nuclear 
fragmentation crossections (o-p, as, Ops) used in the present work are taken from Webber et al. 1992 & 1998 and are listed in Table 2. We 
assume the ISM target density which consists mainly of hydrogen atom to be n = 1 cm^'"* and we take into account the solar modulation 
effect by taking the modulation parameter to be $ = 500 MV. Note that the necessity for the break in the diffusion coefficient at Eo to 
reproduce the peak in the observed data somewhere around 1 GeV/n is consistent with the earlier studies based on the diffusive and leaky 
box propagation models. It is also worth mentioning that in the case of reacceleration models, the peak can be explained using a single 
power-law diffusion coefficient without assuming any break in energy (see e.g., Seo & Ptuskin 1994). However, the observed decrease of the 
secondary abundances with energy above about 1 GeV/n suggests that reacceleration in the ISM cannot be regarded as the dominant process 
for particles with energies > 1 GeV/n (see e.g., Hayakawa 1969). Therefore, by neglecting particle reacceleration as well as other possible 
low energy effects like convection and energy losses, we assume that Eq. (1) properly describes the propagation of CRs in the Galaxy at least 
for energies greater than ~ 1 GeV/n. 

Figure 2 (left) shows the expected variations of ^^C densities at the Earth due to the presence of a nearby discrete source. We choose 
the distance to the source as 0.2 kpc and the source spectral index as P = 2.25. The density variations are obtained from the ratio N'''" /N''^ 
where N''" and A'^*'* denote the densities due to the discrete (given by Eq. 2) and the background (Eq. 12) sources respectively. In Figure 2 
(left), the curves from right to left represent for energies (10 — 10^) GeV/n. We can see that higher energy particles can produce variations of 
much larger amplitudes than the lower energy ones. This can be understood from Figure 2 (right) where we have plotted only the maximum 
density variations (represented by the solid squares) at different energies, i.e. the peak values for each of the curves shown in Figure 2 (left). 
Neglecting particle losses due to nuclear fragmentations, the maximum density of particles of energy E due to a point source located at a 
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Figure 2. Left: ^ C temporal density variations expected at the Earth for energies (10 — 10') GeV/n due to the presence of a nearby discrete source. For the 
calculation, we assume the source distance r = 0.2 kpc, to = 0, F = 2.25, 3? = 25 Myr^^ kpc^^ and <J> = 500 MV. The variations are calculated from the 
ratio N'^'' /N^^ , where A'^'** and N^^ denote the CR densities due to the discrete and the background sources respectively. Right: Maximum density variations 
(solid squares) expected at different energies, i.e. the peak values for the curves shown in the left panel. The dashed straight line in the figure corresponds 
to B"'' which shows that the maximum density variation increases with energy with a slope close to the diffusion coefficient index. f\f'is{max) denotes the 
maximum density the source can produce at the Earth (see text for details). 
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Figure 3. Left: Calculated C spectra normalized to the observed data at 10 GeV/n. Right: O spectra calculated for the O/C source abundance ratio of 1.4. 
Other model parameters: to = 0, F = 2.25, Si = 25 Myr^^ kpc^^, <I> = 500 MV. The dotted lines represent the background CRs and the solid lines the 
total flux which also include the contributions of the nearby SNRs listed in Table 1 . Data points are taken from the results of different experiments given in Zei 
et al. 2007. 

distance r can be found at an age tmax{E) = r'"'/6Dp{E) and is given by Af''^'""'^^(_B) ex Qp(E). The background CRs as given by Eq. 
(12) follows N''^{E) oc Qp{E)/Dp{E). Therefore, the maximum deviation that a discrete source can produce can be obtained as 

T\rds(max) 

^V5^ « DAE) (16) 

For reference, we have also drawn a straight line (ex S"'') in Figure 2 (right) which shows that the maximum density variation increases 
with energy with a slope equal to the index of the CR diffusion coefficient. The deviation from the straight line at lower energies is due to 
the increasing importance of the nuclear fragmentations compared to the diffusion processes at these energies. Therefore, if a nearby source 
is present in the local ISM, we should expect larger density variations for particles which diffuse faster in the Galaxy. However in actual 
practice, since the observations are made at a particular time and location, the actual variation at an energy E depends on the age(s) and 
distance(s) of the nearby source(s). 



6 APPLICATION TO THE NEARBY SNRS 

In this section, we will investigate the possible effects of the presence of nearby SNRs on the observed primary and secondary spectra and 
also, the subsequent effects on the s/p ratios. For the s/p ratio study, we choose the B/C ratio as an example since it is found to be the most 
well-measured ratio among all the available s/p ratios. 

In our study, we consider only the ^^C and ^®0 primaries since boron secondaries are found to be predominantly produced by the nuclear 
interactions of these two species with the ISM. They contribute roughly 50% and 25% respectively of the overall boron produced (Webber et 
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Figure 4. ^^B secondary spectra at the Earth calculated using the ^^C and ^''O spectra shown in Figure 3. The background flux (dotted line) almost overlap 
with the total flux (solid line) since the total contribution of the nearby SNRs is almost neglible. Solar modulation parameter <I> = 500 MV. Data points: 
Swordy et al. 1990 and Engelmann et al. 1990. 
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Figure 5. Same as the ^•^C spectra shown in Figure 3 (left), but now showing the contributions from the individual SNRs. Dotted line: Background spectra. 
Thick solid line: Total spectra. The thin solid lines labelled as (1 — 13) represent the contributions from the individual SNRs listed in Table 1: 1— SN185, 2— 
G73.9+0.9, 3- HB9, 4- S147, 5- CTAl, 6- G65.3+5.7, 7- G299.2-2.9, 8- HB21, 9- Gl 14.3+0.3, 10- Cygnus Loop, 11- Vela, 12- Monogem and 
13— Geminga. 



al. 1998). The calculated spectra (total as well as background) for the CR primaries and secondaries are shown in Figures 3&:4 respectively 
along with the experimental data. The total spectra represented by the solid lines in the figures correspond to the total CR densities which are 
given by 

Nc =N''J + Y^ N^f- (17) 



No = A^o + I] N[ 



Ne 



K% + Nh% + E ^CB + E <B 



jbg 



jbg 



(18) 
(19) 



Table 2. Fragmentation crossections used in our calculations (taken from Webber et al. 1992 & 1998). The subscripts C',0 &c B denote the ^^C, ^®0 & ^^B 
nuclei respectively and CB & OB represent the ^^C^^^B & ^''O^^^B processes respectively. 



Crossection Value (mbam) 
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(TO 308 

<TB 232 

acB 76.8 

aoB 38.5 
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Figure 6. Expected B/C ratios at the Earth in the presence of the nearby SNRs listed in Table 1, calculated for different particle injection times to = 
(0 — 5 X 10^) yr as indicated in the key box. Experimental data are the same as in Figure 1. For reference, we have also shown the B/C ratio arising purely 
from the CR background (solid line). Model parameters: 0/C = 1.4, T = 2.25, SR = 25 Myr^^ kpc"^ and <J> = 500 MV. 



where the subscripts C,0 and B denote the ^^C, 



^^'O and ^^B nuclei respectively, the superscripts 



and ds denote the background 



(represented by the dotted lines in the figures) and the discrete components respectively, and the summations are over the discrete sources 
i listed in Table 1. The subscripts CB and OB represent the ^^C^^^B and ^''O^^^B processes respectively. We adopt a particle release 
time of to = for the discrete sources. The source spectral index is taken as F = 2.25 and the value of the source normalization constant k 
is chosen such that the resulting total ^^C spectrum is normalized to the observed spectrum at 10 GeV/n for the supernova explosion rate of 
K = 25 Myr~^ kpc~^ in the Galaxy. This is shown in Figure 3 (left). The ^^0 spectra calculated using the source abundance ratio of O/C 
= 1.4 is shown in Figure 3 (right). Note that Engelmann et al. 1990 had used a source abundance ratio of 1.24 to reproduce the observed 
C and O abundances using the simple leaky box CR propagation model. The experimental data in Figure 3 are taken from the compilation 
of different experiments given in Zei et al. 2007. The secondary boron spectrum, calculated using the ^^C and ^®0 primary spectra shown 
in Figure 3, is found to explain the observed boron data quite well. This is shown in Figure 4 where the experimental data are taken from 
Swordy et al. 1990 and Engelmann et al. 1990. 

For the primaries, one can notice from Figure 3 that the inclusion of the nearby SNRs in the study produces noticeable deviations from 
the background flux for energies above ~ 100 GeV/n. For instance, the deviation in the case of ^ C reaches a value of --^ 45% at energies 
around 10 GeV/n. But, from Figure 4 we can see that the deviation in the case of ^^B which is considered here as a purely secondary particle 
is almost negligible. The results for the secondaries obtained in this paper agree quite well with the earlier findings of Biisching et al. 2005, 
but some major differences can be seen in the case of the primaries. They found a typical primary amplitude variations of ^ 20% at almost 
all the energies whereas we find almost no variations upto ~ 100 GeV/n and beyond that we see a slow increase reaching a value of ~ 45% 
at ~ 10"" GeV/n. This can be clearly understood if one examines Figure 5 in detail where we have plotted the individual components of the 
overall ^^C spectrum previously shown in Figure 3 (left). The contributions from the individual SNRs are labelled as 1 — 13 (see the figure 
caption for details). It can be seen that among the SNRs, only the Monogem and the Geminga give the highest contribution below ~ 100 
GeV/n while the rest of the SNRs contribute mostly above > 100 GeV/n. Since the CR flux from a point source depends strongly on the 
age and distance of the source, the low energy CRs from the nearby SNRs have not yet reached us effectively except those coming from 
the Monogem and Geminga SNRs. But still, the maximum contributions of Monogem and Geminga below ~ 100 GeV/n are found to be 
approximately two orders of magnitude less than the overall flux. Above > 100 GeV/n, the major contributors are Vela, G299.2-2.9 and 
SN185, each one of them contributing around 10% of the total flux at different energy regions. Adding the contributions of the other SNRs 
also, the total contribution from the nearby known SNRs comes to ~ 31% of the overall CR flux which corresponds to ~ 45% deviation 
from the background level at energies ~ (lO"* — 10^) GeV/n. 

The deviation in the ^^C primary spectra is expected to give a direct implication on the total B/C ratio. This is shown in Figure 6 where 
we have plotted the ratios (Nb/Nc) expected in the presence of the nearby SNRs. For reference, we have also plotted the ratio expected 
purely from the CR background (solid line). The experimental data are taken from Swordy et al. 1990, Panov et al. 2007 and the compilations 
of different experiments given in Stephens & Streitmatter 1998. Strictly speaking, it is not the age t alone which determines the contribution 
of a discrete source but the propagation time At = t — to of the particles after they are released from the source. But, the value of to is not 
exactly known. It may be even different for different sources though we assume the same value for all the sources in the present work. So, 
in Figure 6 we have considered various particle release times i.e. to = (0, 10'*, 5 x 10"*, 10* & 5 x 10*) yr. Below --^ 100 GeV, it can be 
seen that at all the fgs the ratios show almost zero deviations from the background ratio. Any deviation, if exists, are seen at energies > 100 
GeV. At to < 10* yr, we see deviations of magnitude ~ (16 - 26)% in the energy range of ~ (10* - 10*^) GeV/n. At to = 10* yr, the 
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Figure 7. Same as Figure 5, but calculated for to = 10^ yr. Among all the SNRs listed in Table 1, only Geminga with an estimated age of ^ 3.4 X 10"' 
yr have released CRs in the ISM. Dot-dashed line: Geminga. Solid line: Total flux. Thick-dashed line: Background flux (not clearly visible because it almost 
overlap with the total flux). 



ratio shows the maximum deviation reaching a value of ~ 52% at --^ 6 x 10'* GeV/n and for io > W yr, say at to = 5 x 10'* yr, the effect 
of the nearby SNRs becomes negligible producing almost zero deviation from the background ratio. It is worth mentioning at this point that 
detailed studies based on the diffusive shock acceleration in SNRs suggest that the highest energy particles start leaving the source region 
already at the beginning of the Sedov phase (Berezhko et al. 1996), but the major fraction of the accelerated CRs remain confined for almost 
around 10'' yr for a typical interstellar hydrogen atom density of n = 1 cm^"^ . We will come to this point later again in the next section while 
discussing about the implications of the results obtained in this paper. Just for the sake of completeness, we have plotted the ^^C spectrum 



for io = 10' yr in Figure 7. Here, except Geminga (being an old SNR with age '^ 
released CRs into the local ISM, thereby, leading to a negligible effect at the Earth. 



3.4 X IO'' yr), all the other nearby SNRs have not yet 



7 DISCUSSIONS AND CONCLUSIONS 

We have studied the effect of the presence of nearby SNRs on the CR primary and secondary spectra at the Earth. We see strong variations 
in the primary spectra and almost no variation in the case of the secondaries. The results for the primaries obtained here are quite different 
from those obtained from the simulation studies given in Erlykin & Wolfendale (2001) and Biisching et al. (2005). We see variations mostly 
above ~ 100 GeV/n whereas they showed significant fluctuations at all the energies. In fact, their results represent the density fluctuations 
that one can expect at any arbitrary location in the Galaxy due to the random nature of supernova explosions both in space as well as time. 
But, as far as the position of the Earth is concerned, the actual variations can be determined only when one incorporates the nearby sources 
in the study. We, therefore, include the nearby known SNRs in our analysis and found that the primary variations obtained at the Earth show 
considerable differences from those predicted using the Monte-Carlo simulations. 

Below ~ 100 GeV/n, we have found that the effect of the nearby SNRs on the B/C ratio is negligible (see Figure 6). This implies that 
we can safely rely on the observed ratio to determine the CR diffusion coefficient in the Galaxy particularly below ~ 100 GeV/n. Above 
this energy the observed data will not give a reliable information about the propagation parameter because of the significant contaminations 
of the background CRs by those coming from the nearby SNRs. One should note that the primary CRs observed at the Earth are liberated 
from sources located within a short distance which is of the order of the vertical halo height H (Taillet & Maurin 2003). This can also be 
understood from Fig. 8 where we have plotted the fraction of ^'^C primaries at the Earth originated within a radial distance r for energies 
(10 — 10'') GeV/n. The calculations are performed for _ff = 5 kpc. We can see from the figure that for energies less than 100 GeV/n, 50% 
(70%) of the total CRs are emitted within a distance of ~ 3 kpc (5 kpc). This shows that even though the B/C ratio below ~ 100 GeV/n can 
give reliable informations about the diffusion coefficient, the small Galactic region scanned by the CRs reaching the Earth allows them to 
carry information only for a small fraction of the whole Galaxy (see also Taillet & Maurin 2003). 

We have also studied the effect of different particle release times (to = — 10'' yr) from the SNRs on the B/C ratio. At io < 10'* yr, 
we have found a deviation of ^ (16 — 26)% from the background ratio for energies ~ (10'* — 10®) GeV/n and at io = 10* yr we see the 
maximum deviation reaching a value of ~ 52% at around 6 x 10'* GeV/n. For io > 10* yr, say at io = 5 x 10* yr, the effect of the nearby 
SNRs becomes almost negligible, thereby, making the observed B/C data a reliable quantity for determining the CR propagation parameters 
at all energies. This conclusion is, in fact, supported by the predictions of the diffusive shock acceleration theories in SNRs which show that 
a major fraction of the accelerated CRs remain confined in SNRs for almost around 10' yr (Berezhko et al. 1996). For such a confinement 
time, we see that except Geminga all the other local SNRs listed in Table 1 look quite young and might not have liberated CR particles in the 
local ISM. This can be understood from Figure 7 where we have plotted the '^C spectra for io — 10^ yr. 

Therefore, we can conclude that if one accepts the theoretical view of CRs confinement in SNRs for upto ~ 10'' yr, the effect of the 
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Figure 8. Fraction 
height H = 5 kpc. 



Radial distance rlkpc] 
of C primaries emitted within the radial distance r from the Earth for various energies (10 — 10 ) GeV/n. We assume the vertical halo 



nearby SNRs on the observed CRs is expected to be quite negligible. Then, one can expect the observed s/p ratio to give a good estimation 
of the CR diffusion coefficient at all energies, but only for the small Galactic region traversed by them before reaching the Earth. 



APPENDIX A: SOLUTION OF THE TIME DEPENDENT DIFFUSION EQUATION WITHOUT BOUNDARIES 



In rectangular coordinates, the Green's function G(r, r', t, t') of Eq. (1) satisfies 



d^G d'^G d^G 



^''la^ + 



9y2 



dz 



2hncap5{z)G 



dG 
dt 



-5{x - x')5{y - y')S{z - z')S{t - t') 



Since the paricles are assumed to be liberated at time t — t', Eq. (Al) for t > t' becomes simply 



, a^ a^ d^G 

^ * dx^ dy^ dz^ 



Fig 

2hncapS(z)G - -^ = 



Taking fourier transforms of Eq. (A2) with respect to x and y, we obtain 

d^G 



DpK'G + Dp 



where K^ = k1 + kl and 



G{kx,x ,ky,y ,z,z ,t,t ) 



Now, taking laplace tranform of Eq. (A3) with respect to t, we have 

Dp^ - \DpK^ + 2hncapS{z)] G + 5{z)e 
oz^ 



BG 

2hnc(jp5{z)G - 7^- = 



dx I dyG{x,x ,y,y' ,z,z' ,t,t')e' ^^ "^ 






sG = 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



whereG(A:a;,a:', ky,y' ,z,z' ,s,t') = J„°° Ge^°*df and we have used the condition that at t = t',G{x,x' ,y,y' , z,z' ,t,t') = 5{x — x')S{y — 
y')S{z — z'). Note that in Eq. (A5), we have set the value of z' equal to since we will be assumming that the CR sources are located on the 
Galactic plane itself. Solving Eq. (A5) for the regions above and below the z — plane by using the proper boundary conditions at 2; = ±00, 
we get 



G{z, s, t') = Go(0, s, i')e-l-lV^'+«/Op 



The continuity equation at z = can be obtained by integrating Eq. (A5) over z around as 

+0 

Dr, 



dG 



dz 



2hncapGo + e''===" +*'="'' e"'* 



-0 



Solving for Go from Eqs. (A6 & A7), we get 



Go(0,s,t') = 



ite^jX -\'ikyy st 

2Dp y/K'2 + s/Dp + 2hncap] 



and substituting it back to Eq. (A6), we get 



Giz,s,t') = 



-\z\^/k^+sJd^ X gik^x' +ikyy' ^~st' 



\2Dpy/K-^+s/Dp 



2hnca, 



(A6) 



(A7) 



(A8) 



(A9) 
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Taking an inverse laplace's transform of Eq. (A9), we get (see Abramovitz & Stegun 1964) 

G(fc,, x', ky,y', z, z',t, t') = ^ <^ - be'"'+'' <*-* 'erfc b^/T^ + / } (AlO) 

2^ Dp I \/7r(i — t') V 2\/t — i' 

where g = \z\/ ^ Dp and 6 = IhncGpjI^fDp. Further taking an inverse fourier transform of Eq. (AlO) gives 

r(x'-x)^ + (i,'-y)^ 



G{x,x\y,y\z,z\t,^) = ^— —, - T, -- be''+' <■'-' ^ eric ( b^t^F + J (All) 



Using the Green's function given by Eq. (Al 1), we can easily obtain the CR density due to an arbitrary source Qp(r, _E, t) as 

/oo rt 

dr' dt'G{r,r',t,t')Qp{r',,t') 

For a point source of the form Qp(r', E, t') — Qp{E)S{r')5{t' — ta), we obtain a solution given by 



O (E)e ['^"p'-*-*o)\ f p-gV[4(t-to)l ^ ,^2,, ,^ / , 

Np{r,E,t) = ^^^^^" ' -\ -be'^+'^'-'^^eric (bVt^ro + 

8nD-^'\t-to) [ yJ-Kit - to \ 



2^/r^, 



(A12) 



(A13) 



where in cylindrical coordinates, we can just write x'^ + y^ = r'^. 



APPENDIX B: SOLUTION OF THE STEADY STATE DIFFUSION EQUATION WITH VERTICAL GALACTIC 
BOUNDARIES 

The Green's function of the steady state diffusion equation [Eq. (10)] in rectangular coordinates satisfies 

Dp (^ + + ^) - 2hncapS{z)G = -S{x - x')S{y - y')5{z - z') (Bl) 

Taking fourier's transform with respect to x and y, we get 

- DpK^G + Dp^ - 2hncap5{z)G = -e'''-'''+"'yy' 5(z) (B2) 



/OO /'OO 

dx / dyG{x,x\y,y' ,z, z')e^ 
-OO >/ — OO 



(B3) 



where K = k^ ~{- ky and 

/"OO /"OO 

-*^y ^-' V ' " 

' — OO J — OO 

Also note that we have assigned 2;' = in Eq. (B2). Solving Eq. (B2) for the regions z > and 2 < by using the boundary conditions at 
z = dzff, we get 

G{k^,x\ky,y',z)=Go{k:.,x',ky,y',o f''^^\^f^-\'\^^ (B4) 

sinh(Aii) 

The continuity equation at 2 = is obtained by integrating Eq. (B2) over z around as 



dG 



dz 



+0 



2hncapGo + e'''-=''+'''y«' = (B5) 



Solving for Go from Eqs. (B4 & B5), we get 



Go{k.,x',ky,y',0) = [^j^^KcotHKH)l2hncap] ^^^^ 



and substituting it back to Eq. (B4), we obtain 



ik^x' + ikyy' <i\-nh\K( n — \z\\] 

0(/c., a; ,liy,y,z)- [275^;^coth(i^i/) + 2hncap] smh{KH) ^ ' 

Taking an inverse fourier's transform of Eq. (B7) gives 

-, ^00 ^00 ik^(x-x')-iky(y-y') ^mh\K(H—\z\)] 

G{x,x',y,y',z) = -^ dkx dky—^-———— X ^^^nm^Jii (B8) 

V , ,y,y, ) 47r2 y_^ J-^ ^ [2DpKcoth{KH) + 2hncap] smh{KH) 

Eq. (B8) is a two dimensional fourier transform which can be easily simplified in the form of Hankel transform by changing the variables as 

kx = Kcos(j> ; ky = Ksincj) ; x — x' = {r — r')cos6 ; y — y' — [r — r')sin^ (B9) 

Then, 

dry, z) = -\-r ^f[K(H-\z\)] ^ j^ _ ^^^ ^g^^^ 



p A) sinh(_R'7y) \Kcoth{KH) 
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where Jo is the Bessel function of order 0. Having calculated the Green's function, the CR density at r = due to a uniform source 
distribution Sp{r' , E) — ^Qp{E)5{z') in the Galactic plane with radial distances between ri and r2 from the Galactic center is given by 



Np{z,E) = 27rK f \'dr'G{r = 0,r',z)Qp{E) 
J n 

This gives, 



(BID 



NA^,E)^^%i^r ^^"^[^^^-l-l^. , .[r.HKr.)-r.HKn)UK (B12) 



2n 

^ " '■' ^ii 1 1 1 1 L r\ r r I 



p JO smh(KH)\Kcoth{KH) + :^^^ 



where Ji represents the Bessel function of order 1 and we have used the standard relation 

1 
K 



'"^' r'MKr')dr' = ^ [r23iiKr2) - r,Ji{Kri)] (B13) 

By setting ri = and r2 = R, Eq. (B12) can be used to find the CR primary density due to all the sources extended upto a radial distance 
R from the Earth. 
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